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It is shown that space-based gravitational wave detectors such as DECIGO and/or Big Bang 
Observer (BBO) will provide us with invaluable information on the cosmic thermal history after 
£SJ ■ inflation and they will be able to determine the reheat temperature Tr provided that it lies in the 

— . ' range preferred by the cosmological gravitino problem, Tr ~ 10 5-9 GeV. Therefore it is strongly 

desired that they will be put into practice as soon as possible. 

Although we can probe physics during inflation in the early Universe hj rather precisely now using observations of 
t-H , the anisotropy in the cosmic microwave background (CMB) radiation [2j, cosmic evolution from the end of inflation 
to the beginning of the primordial big-bang nucleosynthesis (BBN) is still in the dark age, when many important 
phenomena such as (p)reheating, baryogenesis, generation and freeze-out of cold dark matter particles, etc. have taken 
O | place. It is desirable to clarify the cosmic thermal history in this regime observationally. 

Here we argue that future space-based gravitational wave detectors such as DECIGO Q or the Big Bang Observer 
qj> (BBO), are very useful for this purpose and may be able to determine the reheat temperature after inflation by 
observing stochastic gravitational radiation background generated during inflation in the frequency range around 
0.1 — 10 Hz where foregrounds from astrophysical objects are separable [4[ 1 - 

We introduce tensor perturbations, hij, around a spatially flat Robertson- Walker metric as 
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ds 2 = -dt 2 + a 2 (t) (Sij + 2h ij )dx i dx j , (1) 
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, with a(t) being the scale factor. Decomposing the tensor metric perturbation to Fourier modes as 
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we find that the two independent degrees of freedom tp behave as two massless minimally coupled scalar fields, where 
efj represents polarization tensor with e^e y = S AA for A, A! = +, x. Applying quantum field theory of a massless 
minimally coupled field in de Sitter spacetime, we find that the Fourier modes are characterized by the following 
vacuum correlation, 

^(t)^(t)) = ^S 3 (k-k')S AA ', (3) 
so that the amplitude per logarithmic frequency interval is given by 

h%(f) = 2(h ij 0{f)) =4x8nG (^M) (4) 



Recently, Population III stars are proposed as a dominant component around the deci-Heltz band as a result of the gravitational 
radiation associated with neutrino emissions However, the amplitude is very uncertain due to the lack of understating on early star 
formation history. This signal would be separable if we adopt reliable abundance of Population III stars and take the duty cycle and 
their angular distribution into account. 
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in the long wavelength regime. Here Mq = (87rG) -1 / 2 is the reduced Planck mass and V[0(/)] is the value of potential 
energy density of the inflaton when the scale corresponding to the frequency / today left the Hubble radius during 
inflation. is an expression for the amplitude of tensor perturbation used frequently in the analysis of cosmological 
observations. 

The above expression gives the initial condition to the solution of each Fourier mode in the post inflationary universe 
which behaves as 

h(f,a) oc a(t)^ ^ P _ pa{t * H(t) ) , * = 2^/00, (5) 

in a power-law background a(t) cx t p with p < 1. Here J n (x) is a Bessel function and ao denotes the current scale 
factor. 

Thus the amplitude of gravitational wave takes a constant value, h(f, a) = until its wave length falls shorter 

than the Hubble radius ff _1 at a — 2irfcio/H = a in (f). From the asymptotic expansion of ([5]), one can see that 
the energy density stored in the tensor perturbation starts to decrease just as radiation once the wavelength falls to 
a sub-horizon scale. Therefore in this regime its relative energy density to the background density remains constant 
during radiation domination, while it decreases when the Universe is dominated by other form of energy such as 
nonrelativistic matter or coherent field oscillation. If such a stage lasts long, energy density in sub-horizon tensor 
perturbations tends to be suppressed and the resultant spectrum is modified from a nearly scale-invariant one. The 
present density parameter of the gravitational radiation per logarithmic frequency interval is described as 

(2tt/) 2 a fa m (f) 



^Gw{fM) = \7r^TK{f)\^r L± ) ■ (6) 



It behaves as ^Gw(f,to) <x /~ 2 (/°) f° r the mode which enters the horizon in the matter (radiation) dominated 
regime. Thus the stochastic gravitational wave background not only carries information on the inflationary regime, 
during which they are generated, but also serves as a probe of the equation of state in the early universe 0] (see also 
[Io| for other applications of stochastic gravitational wave background ). 

Besides direct observation by space laser interferometers, cosmological stochastic gravitational wave background 
can be detected indirectly through CMB. By measuring its B-mode polarization, we can probe the amplitude of tensor 
perturbation from inflation. 

In a simple single-field slow-roll inflation model with a potential V[<j>], observable quantities such as the amplitude of 
curvature perturbation, A-jj, its spectral index, n s , and running, dn s /d\nk, are described by the slow-roll parameters, 

as 
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= 24^%' n > = 1 - 6£ + 2 * ShS = 16£ " - 246 - 

where each quantity should be evaluated at the time of horizon crossing during inflation. 

WMAP has measured both An and n s with an unprecedented accuracy. However, since An depends not only on 
V[<j>] but also on the slow-roll parameter e, we have not been able to fix the energy scale of inflation yet. Detection 
of tensor perturbation is essentially important to determine the energy scale of inflation. 

Theoretically, different inflation models predict different energy scale of accelerated expansion. Among them, chaotic 
inflation model [ll| with a massive scalar potential, VW[ = m 2 (j> 2 /2, which is attractive both from naturalness of 
initial condition and phenomenological point of view [12j |. we find a relatively large amplitude of tensor perturbation: 
r = A 2 t /A 2 J (fc ) = 0.16, where pivot scale, fc , corresponds to the CMB scale. Other models such as small field models 
may predict much smaller value of r including a vanishingly small one, which would make us desperate with regards 
to the detection of B-mode polarization. 

Recently, however, it has been claimed that, if we use the observed value of the scalar spectral index n s = 0.961 ± 
0.017 by WMAP as a constraint, even small-field models would predict r > 10 -3 [HI]. That is, it would be unusual 
to have r/ ^> e because it would mean the inflaton is near an inflection point when the CMB scale left the Hubble 
radius and it would be very difficult to achieve the right number of e-folds thereafter without severe fine tuning. Then 
the observed deviation of the spectral index from unity, n s ~ 0.95 — 0.99, implies e = O(0. 01) or r ~ 0.1. Indeed 
authors of [l3| have considered a number of potentials with different parameters and calculated n s and r under the 
condition that the duration of inflation takes a proper value. As a result they find that the observed value of n s 
implies r > 0.003. 



3 



The same issue has been studied by Boyle, Steinhardt, and Turok [1J| in somewhat more model-independent 
manner. They count the number of zeros of the slow roll parameters e and r\ in the last 60 e-folds of inflation as a 
conservative measure of how many derivatives of them must be fine tuned to achieve a given set of (n s , r) . They find 
inflation models with no fine tuning give n s < 0.98 and r > 10~ 2 and that models predicting r < 10 -3 require nine 
or more extra degrees of fine tuning. 

If r is indeed larger than, say, 0.003, we will be able to detect B-mode polarization by ongoing and planned projects 
and its implication is profound. Because we already know that ~ 2.0 x 10~ 9 on CMB scale @, by measuring r 
we can fix the energy scale of inflation as 

V[<t>] = (3.2 x 10 16 GeV) 4 r = (7.5 x 10 15 GeV) 4 (^) . 

Once we succeed in measuring r by B-mode polarization, we will be even more confident in simple slow-roll single- 
field inflation, on which the above arguments jl3l . [Tij are based, and we can predict the amplitude of gravitational 
waves on frequencies accessible by the other means of observation, namely, a space laser interferometer such as 
DECIGO [H and BBO. Extrapolating the amplitude of tensor perturbation on the CMB scale with the wavenumber 
k = 0.002Mpc~ 1 corresponding to the frequency //, = 3 X 10 _18 Hz today, to higher frequency using the slow-roll 
parameters measured at fco we find 



A 2 h (f)=Al(f h ) 



2eln ±- + 2e(r 7 - e) In ±-\ + -e(-12e" + 16ery - ±rf - 2£) In f- 
Jh \ Jh I 3 V Jh 



(7) 



Among the planned space laser interferometers, LISA [15j is most sensitive to frequency around 10 -3 Hz where 
stochastic gravitational wave background is dominated by astronomical sources such as white dwarf binaries and 
it seems difficult to detect inflationary gravitational wave background even if r is maximal. On the other hand, 
DECIGO or BBO targets frequency window around / = 0.1 — 1 Hz where contamination of astrophysical foregrounds 
is separable, so they are ideal to probe inflation. In the standard cosmology the gravitational wave with its current 
frequency 0.1 Hz reentered the Hubble radius in the radiation dominated regime when the cosmic temperature was 
T = 4 x 10 6 GeV assuming that there were no significant entropy production after that. 

Interestingly, this temperature is in the expected range of the upper bound on the reheat temperature, Tr, after 
inflation imposed by the decay of unstable gravitinos which are produced by thermal scattering in the reheating 
stage. Specifically, in order to ensure successful BBN in the presence of the hadronic decay of gravitinos, the reheat 
temperature should satisfy Tr < 10 6_8 GeV depending on the gravitino mass and the hadronic branching ratio [l7l | . 
In case gravitino is the lightest supersymmetric particle and hence stable due to R-parity conservation, we also find 
a constraint on the reheat temperature so that it does not overdose the universe, 

for the gravitino mass m 3 / 2 = 10~ 4 — 10 GeV [Hj], where rrig denotes the gluino mass. 

If the reheat temperature indeed satisfies the above constraints, it is likely that the universe is still dominated by 
coherent inflaton field oscillation when the frequency range probed by DECIGO/BBO reentered the Hubble radius. 
Since the cosmic expansion law is the same as in the matter-dominated era as field oscillation is driven by a mass 
term, the resultant spectrum of £Igw{J , io) acquires a modulation proportional to f~ 2 for those frequency bands 
entering the horizon in this regime. 

So far we have implicitly assumed that there is no significant entropy production after the completion of reheating 
after inflation. If nonncgligible amount of entropy is produced in the lower temperature regime, from the decay of a 
long-lived scalar field other than the inflaton for example, not only the spectrum of QgwU >^o) but also constraint 
on the reheat temperature imposed by the thermal gravitino problem are modified. Defining the entropy increase 
factor F by the ratio of entropy in a fixed comoving volume before and after the entropy production Q , we find that 
the overall amplitude of the high-frequency part of Qowifito), the frequency of gravitational radiation entering the 
horizon at the end of reheating, and the constraint on the reheat temperature by the gravitino problem are multiplied 
by F -4 / 3 , _F -1 / 3 , and F, respectively. 

Taking all these factors into account we can numerically calculate the spectrum of stochastic gravitational wave 
background in the band probed by DECIGO/BBO using future-observed values of the tensor-to-scalar ratio and the 
slow-roll parameters. Figure 1 summarizes the result of numerical calculation for the case r = 0.16, e = 77 = 0.01, and 
£ = 0, which are realized in chaotic inflation driven by a quadratic potential. 
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FIG. 1: Density parameter of gravitational radiation for different reheat temperature Tr and entropy increase factor F in 
chaotic inflation with a massive scalar field. Curves without F imply F = 1. The region above the thick curve with / > 0.1Hz 
can be observable by the ultimate DECIGO after its ten years' operation. 



We find a nearly flat spectrum 



ficw(/,*o) = 2.8 x 10- 16 F- 4 /-Vio6 3 7 5 (^y) 
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in the low frequency band reentering the horizon in the radiation dominated regime after reheating (and before 
possible entropy production), 



/ < 0.042 



TrF- 1 / 3 
10 7 GeV 



Hz = /_ 



(10) 



Here 3*106.75 denotes effective number of relativistic degrees of freedom normalized by 106.75 when the relevant mode 
reentered the Hubble radius. For higher frequency modes with 



we obtain a nearly power-law spectrum 



h 



0.31 



h_ 

f 

10 7 GeV 



Hz. 



(11) 

(12) 
(13) 



The cosmological information we can obtain by observing these gravitational wave background depends on the shape 
of the spectrum observable at the DECIGO/BBO band, namely, 0.1Hz < / < 10 Hz, which can be classified to the 
following three cases (a)-(c). 
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Case (a) /_ > 10 Hz corresponding to Tr > 2.4 x lO 9 ^ 1 / 3 GeV: We can observe the nearly-flat region of the 
spectrum and determine F from the overall amplitude of £Igw- Since WMAP constrains r < 0.55 0, F should satisfy 
F < 3 x 10 4 . Thus for stable gravitino, the above lower bound on Tr indicates 

and for unstable one, larger mass (ro 3 / 2 > 10 TeV) is preferred. 

Case (b) /_ < 10 Hz and /+ > 0.1 Hz corresponding to 7.1 x lO 5 ^ 1 / 3 GeV < T R < 2.4 x lO 9 ^ 1 / 3 GeV: We can 
fix F from the overall amplitude and Tr from the shape of the spectrum independently. This is the ideal case that 
space laser interferometers can practically determine the entire thermal history of the universe between inflation and 
BBN. 

Case (c) /_|_ < 0.1 Hz corresponding to Tr < 7.1 x 10 5 F 1 / 3 GeV: We can observe only the power-law region of the 
spectrum and measure only the ratio Tr/F. This ratio, however, fixes the gravitino-to-entropy ratio uniquely apart 
from a logarithmic correction. 

In conclusion, space laser interferometers such as DECIGO and BBO will bring about invaluable information on 
the delayed reheating stage required by the gravitino problem after inflation and will be able to determine the reheat 
temperature and/or the entropy increase factor. Hence it is desired that they are put into practice as soon as possible. 
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